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Abstract 



We study resonant tunnelling through double-barrier structures under an 
applied bias voltage, in which nonlinearities due to self-interaction of electrons 
in the barrier regions are included. As an approximation, we concern ourselves 
with thin barriers simulated by 5-function potentials. This approximation 
allows for an analytical expression of the transmission probability through 
the structure. We show that the typical peaks due to resonant tunneling 
decrease and broaden as nonlinearity increases. The main conclusion is that 
nonlinear effects degrade the peak-to-valley ratio but improve the maximum 
operation frequency of the resonant tunnelling devices. 
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I. INTRODUCTION 



Resonant tunnelling (RT) through a heterostructure quantum well has recently attracted 
considerable attention because of its applications as ultra-high-speed electronic devices. As 
an example, RT in GaAs/GaAlAs double-barriers at THz frecuencies has been reported in 
the literature [I]. The way RT arises is a basic phenomenom which can be explained with 
elementary quantum mechanics [0]: there exists a dramatic increase of the transmission 
probability whenever the energy of an incident electron matches one of the unoccupied quasi- 
bound states inside the well. The energy of the incident electron is close to the Fermi energy, 
which is constant in the device. Therefore, to obtain a RT current in practice, one must 
modify the energy of the discrete level in the well applying a bias voltage. If initially there 
are no electrons in the well and the electrons are forced to tunnel by the applied voltage, after 
a time of a few ~ H/T, T being the width of the discrete level, high transmission is achieved 
by multiple reflection. This mechanism is physically analogous to what occurs in a Fabry- 
Perot resonator. Since multiple scattering of waves is required for RT, it can be realized that 
such a phenomenon depends crucially on the linear superposition principle. Then a natural 
question arises in this context, namely how spatial nonlinearities affect RT. Nonlinearity 
may arise in semiconductors when electrons polarize the surrounding medium, which reacts 
to change the electron state by a fedback process, or when electron-electron interaction 
is taken into account. There are two important magnitudes which, in principle, could be 
largely modified by nonlinearity: the peak-to-valley ratio and the level width T, which is 
related to the reliability of the device at high frequencies. For instance, several satellite 
peaks in addition to the main peak have recently been found via coupled electron-phonon 
states in RT double-barrier heterostructures ||. Moreover, the occurrence of nonlinearities 
opens the possibility of studying new and interesting phenomena in semiconductor devices, 
like bist ability or chaotic wavef unctions [[| 

In this Letter we focus our attention on a simple model to study nonlinear effects in 
double-barrier heterostructures. We consider very narrow barriers of a nonlinear material 
embedded in a linear material. We simulate the two barriers by means of 5-function poten- 
tials, which allow for an analytical treatment even in the presence of nonlinearities. Similar 
models have been proposed to study wave propagation in nonlinear periodic media f| and 
nonlinear Stark- Wannier resonances . Thus, the propagation of a stationary wave in the 
nonlinear double-barrier heterostructure is described by means of the following Schrodinger 
equation (we use units such that H — m = 1) 

- tf"(a;) + [5(x) + 5(x - a)} [(3 + a*(x)**(x)] + V(x)V(x) = E (1) 

The parameter {3 is the strength of the (^-potential in the linear case (i. e. the area of the 
semiconductor barrier) and a measures the strength of the nonlinear coupling. The width 
of the quantum well is denoted by a. The potential due to the applied bias is V(x). This 
potential may be simulated by step functions, which have been demonstrated to be a good 
approach to linearly rising potentials, reproducing accurately the properties of the electronic 
states, even the subtle ones J7J. Therefore we take 

f 0, x<0, 

V(x) = | - V, 0<x<a } (2) 
' —V, x > a, 
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where V is the applied voltage. 

The solution of Eq. (fj) in the linear semiconductor is a combination of travelling waves. 
As usual in scattering problems, we assume an electron incident from the left and evaluate 
the transmission and reflection probabilities when passing through the structure. To this 
end we take 

A (e ikQX + re~ ikox ) , x < 0, 
*(x) = { A [Ce iklX + De~ iklX ) < x < a, (3) 
Ate ik2X x > a. 

t and r denote the transmission and reflection amplitudes, and we have defined for simplicity 
fcg = E, k\ = E + V/2, and k 2 — E + V. The transmission probability is then found to be 
T = (ko/k2)\t\ 2 . To compute it we must set appropriate boundary conditions at the barriers, 
in order to eliminate the two unknowns C and D. These can be obtained by integrating 
Eq. (|1]) around x = and x = a. In so doing, we find that the wave function itself is 
continuous whereas its derivative satisfies 

tf'(0 + ) - tf'((T) = [/? + a|^(0)| 2 l tf(0), 



f3 + a\^(a)\ 2 *(a) (4) 



Using these boundary conditions and solution @ we have after a little algebra 

r = t jcos£;ia+ — ^ — [(f3 + a\t\ 2 ) - ik 2 ]\ - 1, (5a) 



sin k\a 

t = (l + r ) Smkia \ \ki cot ha + ik + (3 + all + r\ 2 } - 2ik r\ . (5b) 
k\ L L J ' 

where a = <5|y4| 2 . Finally, inserting (pa[) in (|5b|) we obtain an expresion for the transmission 
amplitude, and thus the transmission probability can be found for a given set of parameters 
(E, a, P, a and V). 

Figure [I] shows the transmission probability as a function of the applied voltage V and the 
nonlinear coupling a. The energy of the particle is taken to be E = 2, (3 = 10, and a = 1.5. It 
is important to mention that the qualitative aspects of the results remain almost unchanged 
when varying the parameters E, (3, and a. From this plot it is clear that several peaks appear 
on incresing the applied voltage, no matter the value of the nonlinear coupling a. Each peak 
corresponds to a quasi-bound state level in the well. It is also apparent that the width of these 
peaks is relatively large, specially for higher order ones. At this point we should remark that 
the width of the peaks appearing in Fig. [1] (transmission probability versus applied voltage) 
is not exactly equal to T, i. e., the level width. Nevertheless, assuming that the lifetime is 
not strongly dependent on the position of the level (so that it remains almost unchanged 
under minor variations of the applied voltage), it is clear that the width of the resonance 
peak in the transmission coefficient is very close to T. Hence we can accurately determine T 
from Fig. 0. The fact that T is relatively large is easily understood by considering that the 
barriers are, in fact, very narrow, so that the coupling of the quantum well states with the 
continuum is very strong. Thus the lifetime is very short and, accordingly, V is large. Also 
r increases upon increasing the order of the quasi-bound state in the well, as expected. Now 
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let us comment the main effects of nonlinearity. In Fig. p] one can see that the peak-to- valley 
ratio decreases dramatically on increasing a from —30 up to 30. Moreover, it is also clear 
that the resonances shift to higher voltages when increasing a > but, on the contrary, 
resonances shift to lower voltages for a < 0. This means that the quasi-bound states within 
the well region are raised on incresing nonlinearity. With more elaborated models, it has 
been recently found a similar shift to higher voltages whenever Coulomb repulsion between 
electrons increases in quasi-one-dimensional systems M. Hence we are led to the conclusion 
that our model retains most of the physics of resonant tunneling, since Coulomb repulsion 
can be regarded as a possible origin of nonlinearities within the one-particle approach. 

Since the width of the peaks is a very important parameter for the design of high-speed 
devices, we have also studied T as a function of the nonlinear parameter. Results are shown 
in Fig. [2] for three different values of the nonlinear parameter a. Besides the relative shift of 
the peaks above mentioned, it is seen that T decreases with a for peaks of the same order. 
Hence we are led to the conclusion that, in order to improve the quality of the RT device, 
one should look for a compromise between the peak-to-valley ratio which decreases with 
a, and the width of the levels which also decrease. Hence, on increasing nonlinearity, one 
obtains a lower RT current but, in principle, a higher operation frequency. 

A possible extension of this work concerns narrow-gap semiconductors for RT P,|TT|. In 
this case RT takes place by interband tunnelling mechanism, and a simple one-band descrip- 
tion is not longer valid. Instead, one must use more elaborated band structure including 
the effects of several bands of the host semiconductors. One of the simplest model is the 



two-band model ||11|| , in which envelope-functions of the conduction- and valence-bands are 
included. The equation of motion is then analogous to the Dirac equation for relativistic 
electrons. Hence, Eq. [I] should be replaced by a Dirac-like equation with nonlinear terms. 
Fortunately, the nonlinear Dirac equation with 5-function potentials is also exactly solvable 



14 and, therefore, the generalization of the present work to narrow-gap heterostructures 



seems to be achievable. In addition, more realistic and elaborated models are indeed required 
to describe nonlinear electron dynamics through thick barriers even if only one envelope func- 
tion is used, for which the 5-function approach is no longer valid. Work in that direction is 
in progress. 
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FIGURES 



FIG. 1. Transmission probability T as a function of the applied voltage V and the nonlinear 
coupling a for E = 2, [3 = 10, and a = 1.5. 

FIG. 2. Transmission probability T as a function of the applied voltage V for E = 2, f3 = 10, 
a = 1.5, and (a) a = -10, (b) 0, and (c) +10. 
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